Finite Volume Analysis of Nonlinear 
Thermo- mechanical Dynamics of 
Shape Memory Alloys 

L. X. Wangi and Roderick V.N. Melnik^* 
^MCI, Faculty of Science and Engineering, 
University of Southern Denmark, 
^ ■ Sonderborg, DK-6400, Denmark 

O ; ^Mathematical Modelling and Computational Sciences, 

; Wilfrid Laurier University, 75 University Ave West, 

: Waterloo, ON, Canada N2L 3C5 

oo 
(N 

pq . Abstract 

In this paper, the finite volume method is developed to analyze coupled dynamic 
^ ' problems of nonlinear thermoelasticity. The major focus is given to the description of 

martensitic phase transformations essential in the modelling of shape memory alloys. 
Computational experiments are carried out to study the thermo-mechanical wave inter- 
^ ■ actions in a shape memory alloy rod, and a patch. Both mechanically and thermally 

. induced phase transformations, as well as hysteresis effects, in a one-dimensional struc- 

I ture are successfully simulated with the developed methodology. In the two-dimensional 

^SJ ' case, the main focus is given to square-to-rectangular transformations and examples of 

. martensitic combinations under different mechanical loadings are provided. 

o 
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^ ' finite volume method. 



X ■ 1 Introduction 



The existing and potential applications of Shape Memory Alloys (SMA) lead to an increasing 
interest to the analysis of these materials by means of both experimental and theoretical ap- 
proaches These materials have unique properties thanks to their unique ability to undergo 
reversible phase transformations when subjected to appropriate thermal and/or mechanical 
loadings. Mathematical modelling tools play an important role in studying such transforma- 
tions and computational experiments, based on mathematical models, can be carried out to 
predict the response of the material under various loadings, different types of phase transfor- 
mations, and reorientations. The development of such tools is far from straightforward even 
in the one-dimensional case where the analysis of the dynamics is quite involved due to a 
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strongly nonlinear pattern of interactions between mechanical and thermal fields (e.g., [HI ITT] 
and references therein). For a number of practical applications a better understanding of the 
dynamics of SMA structures with dimensions higher than one becomes critical. This makes the 
investigation more demanding, both theoretically and numerically. 

Most results reported so far for the one-dimensional case have been obtained with the Finite 
Element Method (FEM) [11EII23]- In addition to the challenges pertinent to coupling effects, 
we have to deal also with strong nonlinearities of the problem at hand. One of the approaches 
is to employ a FEM using cubic spline basis functions, in which case the nonlinear terms can 
be smoothed out by one of the available averaging algorithms. As an explicit time integration 
is typically employed in such situations, this results in a very small time step discretization. 
Seeking for a more efficient numerical approach, Melnik et al. fiSi [20] used a differential- 
algebraic methodology to study the dynamics of martensitic transformations in a SMA rod. 
An extension of that approach has been recently developed in [121 dZl where the authors 
constructed a fully conservative, second-order finite-difference scheme that allowed them to 
carry out computations on a minimal stencil. However, a direct generalization of the scheme 
to a higher dimensional case appeared to be difficult. 

In this paper, we approach the same problem from the Finite Volume Method (FVM) point 
of view. The method is based on the integral form of the governing equations, leading to 
inherently conservative properties of FVM numerical schemes. The methodology is well suited 
for treating complicated, coupled multiphysics nonlinear problems [H [6], [7] . It can be relatively 
easily generalized to higher dimensional cases. In addition to its wide-spread popularity in 
CFD, the method has been applied previously to linear elastic and thermoelastic problems 
[21 ini [3, [IS]- There are several recent results on the application of FVM to nonlinear thermo- 
mechanical problems and nonlinear elastic problems [11 [25]. In this paper, we develop a FVM 
specifically in the context of studying martensitic transformations in SMAs and demonstrate 
its performance in simulating the dynamical behavior of SMA rods and patches. 

The paper is organized as follows. The mathematical models for the dynamics of marten- 
sitic transformations in ID and 2D SMA structures are described in Section 2. Key issues 
of numerical discretization of these models, including the FVM and its computational imple- 
mentation via the Differential-Algebraic Equations (DAE) approach, are discussed in Section 
3. Mechanically and thermally induced transformations and hysteresis effects in SMA rods are 
analyzed in Section 4. Section 5 is devoted to studying nonlinear thermomechanical behavior 
and square-to-rectangular transformations in a SMA patch. Finally, conclusions are given in 
Section 6. 

2 Mathematical Model for SMA Dynamics 

We start our consideration from a mathematical model for the SMA dynamics based on a 
coupled system of the three fundamental laws, conservation of mass, linear momentum, and 
energy balance, in a way we described previously in [211 El [26]. Using these laws, the system 
that describes coupled thermo-mechanical wave interactions for the first order martensitic phase 
transformations in a three dimensional SMA structure can be written as follows [TBI [211 El] 




z,j = 1,2,3 



(1) 
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where p is the density of the material, u = {Mi}|i=i,2,3 is the displacement vector, v is the 
velocity, cr = {cxjj} is the stress tensor, q is the heat flux, e is the internal energy, f = (/i, /2, fzY 
and g are distributed mechanical and thermal loadings, respectively. Let cf) be the free energy 
function of a thermo-mechanical system described by ([T]), then, the stress and the internal 
energy function are connected with cj) by the following relationships: 

where 9 is the temperature, and 77 the Cauchy-Lagrangian strain tensor defined as follows: 

In what follows, we employ the Landau- Ginzburg form of the free energy function for both 
ID and 2D SMA dynamical models [H [HI HE]- In the 2D case, we focus our attention on the 
square-to-rectangular transformations that can be regarded as a 2D analog of the realistic cubic- 
to-tetragonal and tetragonal-to-orthorhombic transformations pTl [T2] . It is known that for this 
kind of transformations, the free energy function can be constructed by taking advantage of 
a Landau free energy function F^. In particular, following [HI [121 US] (see also references 
therein), we have: 

(f) = -c^e \n.e + ]:aie\ + ^0363 + F^, Fl = {9 - 9q) el - ^0463 + ^a^el, (4) 
II / 4 

where c„ is the specific heat constant, 9q is the reference temperature for the martensite tran- 
sition, Oj, i = 1,2,3,4,6 are the material-specific coefficients, and ei, 62, 63 are dilatational, 
deviatoric, and shear components of strain, respectively. The latter are defined as follows: 



ei = (r^ii + 7722) /v2, 62 = (r^ii - 7722) /v2, = {tju + V21) f^- (5) 

This free energy function is a convex function of the chosen order parameters when the tem- 
perature is much higher than 9o, in which case only austenite is stable. When the temperature 
is much lower than 9q, becomes non-convex and has two local minima associated with two 
martensite variants, which are the only stable variants. If the temperature is around 9q, the 
free energy function has totally three local minima, two of which are symmetric and associated 
with the martensitic phases and the remaining one is associated with the austenitic phase. In 
this case both martensite and austenite phases could co-exist in the system. By substituting 
the above free energy function into the conservation laws for momentum and energy, and using 
Fourier's heat fiux definition 

q = -k9^ (6) 

with k > being the heat conductivity of the material, the governing equations for 2D SMA 
patches can be written in the following form: 

d^^l V2 d . 3 5\ d fl \ ^. 

P~di^ ^ ~Ydx ^'^^^^ ^ '^^^ ~ '^^^^ ^^^^^ dy [2^^^^ ) 

d^U2 d fl \ V2 d , ,^ ^ , 3 5. ^ 

89 , fd^9 d^9\ ^ de2 



(7) 
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As always, we complete system ([7]) by appropriate initial and boundary conditions which are 
problem specific (see Sections 4 and 5). As discussed before in [iniES], the 2D model given by 
EqJT] can be reduced to the Falk model in the ID case 

ot ox'^ ox ot 

where fci, ^3, Cv and k are re- normalized material-specific constants, 9i is the reference 
temperature for ID martensitic transformations, and F and G are distributed mechanical and 
thermal loadings. 

In the subsequent sections, the above models are applied to the description of the first order 
martensitic transformations. While such transformations are reasonably well documented for 
the ID case, only few results are known for the 2D case. In what follows, we develop a FVM 
to simulate the dynamics described by the models ([7]) and ([H]) and apply it in both ID and 2D 
cases. 



3 Numerical Algorithm 

The systems ([7]) and ([H]) are analyzed numerically with the FVM implemented here with the 
help of the DAE approach. For the ID case, the FVM method yields the same result as 
the conservative scheme already discussed in [17]. However, the approach developed here is 
generalized in a straightforward manner to a higher dimensional case and we demonstrate its 
applicability by a numerical example in the case of two spatial dimensions. First, we note that 
it is convenient to replace the original model ([7]) by a system of equivalent differential-algebraic 
equations as it was proposed earlier in 



dei a/2 / dvi dv2\ de2 \/2 / dvi dv2 
dt 2 \dx dy J ' dt 2 \dx dy 
dvi dan dai2 „ dv2 dau da 



P— = + + P— = + + Jy 

(9) 



dt dx ' dy ' ^ dt dx ' dy ' 
de , fd'^e d^9\ ^ de2 



V2 

an = —{aiei + a2{9 - Oq) 62 - 0462 + a^el), 
1 

(^12 = 0-21 = -^363, 
V2 

(^22 = ^(cnei - a2{6 - 9q) 62 + a^el - a^el). 

This system is solved numerically together with the compatibility relation written below in 
terms of strains: 

d'^ei d'^ei r- d'^e^ 9^62 9^62 



dx\ dx\ ^ dxidx2 dx\ ~'~ dx\ ^ ^ 
There are eight variables in total that the problem needs to be solved for in this 2D case and 
there are eight equations. The equations for strains, velocities and temperature are all differ- 
ential equations, complemented by stress-strain relationships which are treated as algebraic. 
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In what follows, we highlight the key elements of our numerical procedure based on the 
FVM implemented with the help of the DAE approach. First, all equations in the system 
are discretized on a staggered grid represented schematically in Fig 1. Assuming that the entire 
computational domain is a rectangle with an area of x Ly cm^, we define the spatial integer 
grid points {xi,yj) and the spatial flux points {xi,yj) as follows: 



Xi = ih^, ,z = 0, 1,2, ■ ■ ■,M, Xi = {i-^)h^, ,z = l,2, ■■■,M 
yj=jhy^ ,J = 0,1,2, ■■■,A^, y. = [j-^)hy, ,j = 1,2,---,N 



where M and N are the number of grid points such that M x hx = and N x hy = Ly, 
respectively. The (z, j)th control volume for the velocities is [a;j,Xj+i] x [yj,yj+i], as sketched 
by the rectangular tiled mosaic area in Fig.l, including the upper right part overlapped with 
the hatched area. The variables, defined in this control volume, that will be differentiated 
are marked by a top bar, for instance Vi{i,j). The control volume for the strains ei and 62, 
temperature 9, and stresses an, au, and a22 is given by [xi, Xi+i] x [yj, yj+i], represented by the 
rectangular hatched area in Fig. 1. We refer to these variables, defined in this control volume, 
without a top bar, for instance ei(i, j) for ei, etc. 

By integrating all the differential equations over their own control volumes and assuming 
that all the unknowns are linear in each single control volume while being continuous and 
piecewise linear in the entire computational domain, the five partial differential equations are 
reduced to a system of ordinary differential equations. The remaining three algebraic equations 
of the original system are discretized directly on the grid. The result is the following system: 

dei{i,3) IT - r ■\ , r n - r ■ 
— — = {IyDxVi{i,j) + IxDyV2{i,j 

= ilyDxMhj) - I.DyV2it,j))/V2, 

P^^^^%^ = 4^x^11 - 1, J - 1) + IxDya^2i^ - 1, J - 1) + /i, 

P^^^^%^ = W^i2(^ - 1, J - 1) + LDya22{t - 1, J - 1) + f2, 

= -^(«iei(«, j) +a2(^(^, j) - ^o)e2(^, j) - ^giihj) + ^g2{hj)), 
(^i2{i,j) = (X2i(2, j) = ^(0363(^,7)), 

(^22{hj) = ^(aieiiij) - a2{9{i,i) - ^o)e2(i, j) + ^gi{i,j) - ^92(1, j))- 

where and Dy are the discrete difference operators in the x and y directions, respectively, 
while Ix and ly are the discrete interpolation operator in the x and y directions, and A is the 
discrete Laplace operator. For example, for the simplest case of the first order accurate scheme, 
the operators and ly could be written as follows 

DxViiiJ) = [viii,] + 1) -vi{i,i))/hx, IyVi{i,j) = {vi{i,j) + Vi{i + lJ))/2, (13) 



(12) 



5 



with similar representations for the second order accurate schemes. Moving to the time dis- 
cretization procedure, it is convenient to re-write system fll2p in the following vector-matrix 
form: 

A^ + H(t,X,U) = (14) 

with matrix A = diag(ai, 02, a^) having entries "one" for differential and "zero" for algebraic 
equations for stress-strain relationships, and vector-function H defined by the right hand side 
parts of ( fT2l) . This (stiff) system is solved with respect to the vector of unknowns U that 
have 6 x rrix x niy + 2 x {rux -|- 1) x {ruy + 1) components by using the second order backward 
differentiation formula (BDF) [TO] : 

A Qu" - 2U'^-i + ^U"~2^ + AtH (t„, X, U") = (15) 

where n denotes the current time layer. 

This spatio-temporal discretization is applied to the analysis of phase transformations with 
the following modification. In order to improve convergence properties of the scheme, we employ 
a relaxation process connecting two consecutive time layers via a relaxation factor u as follows: 

y{z,j) = (1 - ^) X y{z,jr + ^ X y(z, j)"+\ (16) 

where the variable y could be any of the following: ei(z,j), 62(1, j), Vi{i,j), V2{i,j), or 0{i,j). 
Note that in the general case the relaxation factors need not be the same for all the variables. In 
the present paper, all the numerical results have been obtained using ffTUl) with all the relaxation 
factors set to 0.85. 

We note that nonlinear terms in the model are averaged in the Steklov sense [IT], so that 
for nonlinear function /(e2) (in particular, for e| and 63), averaged in the interval [63 , 62"*"^], we 
have 

1 f 2 



"^2 

3 o -r-i ^5 



9 {el e^i) = / fie2)de2. (17) 



Applying this idea to 63 and e^, we have: 



91 [i,J) = T^T^ = 77 .16; 



^2) _ J- -^-^ ,n+l\3-k f n\k 



k=0 

+l\5—kfn\k 



where 63 and 63^^ stands for e2(i, j)" and 62(1, respectively. 

Finally, we note that in our FVM implementation the nonlinear coupling term in the energy 
balance equation is regarded as a time-dependent source term. In the (i,j)th control volume 
for the discretization of 6, we approximate that term as follows: 

hee2^ dxdy ^ e{z,j)e2{z,j)^^^. (19) 



dt J ^ ' dt 



As seen from ( 1T5|) . we use an implicit time integrator based on the BDF. At each time step 
we apply the bi-conjugate gradient method to solve the resultant system of algebraic equations 
with the Jacobian matrix updated on each iteration. 
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4 Dynamics of SMA Rods and Strips 



We first consider a situation where the deformation of a 2D SMA sample in the xi direction 
substantially exceeds the deformation in the other direction, so that the deformation in the X2 
direction can be neglected and the sample can be treated as a SMA long strip or simply as a 
rod. Introducing formally e = du/dx and v = du/dt, system can be recast in the following 
form: 

de dv dv ds ^ 
dt dx' ^ dt dx 
s = ki{e -ei)e- k2e^ + fc3e^ (20) 

de a^e , 8v ^ 

ot ox^ ox 

where e is strain and s is stress. 

The numerical procedure described in Section 3 is applied here to the solution of system fpUj) . 
It is aimed at the analysis of martensitic transformations in the SMA rod, including hysteresis 
effects during the transformations. Computational experiments reported in this section were 
performed for a Au2zCuzoZn4rj rod with a length of L = 1 cm and all parameter values found 
in [HlIiniES], in particular: 

ki = 480 g/ms^anK, A;2 = 6 x lO^g/ms'^cmK, k^ = 4.5 x 10^ g/ms'^cmK, 

01 = 208K, p = ll.lg/cm^, = 3.127Ag /ms'^cmK, k = 1.9 x 10~^ cmg/ms^K. 

The boundary conditions for u and 6 for all the numerical experiments reported in this 
section are: 

06 do 
u(0,t) = Ui(t), umt) = unit),qquad—{0,t) = eL{t), —{L,t) = eR{t) (21) 

with given functions Uiit) and 6'j(t), i = L,R and corresponding conditions for the velocities. 

In the numerical experiments reported below, we used only 9 nodes for the velocity dis- 
cretization (and 8, excluding boundaries, for the rest of variables). The time stepsize in all the 
experiments was set to r = 1. x 10~^. All the simulations were performed for the time period 
[0, 24] which spans two periods of the loading cycle. 



4.1 Mechanically Induced Transformations and Hysteresis 

The first numerical experiment deals with the case of mechanical loading in the low-temperature 
regime. The initial conditions for this computational experiment are defined by the following 
configuration of martensites ([IH [ISl 121]): 

( 0.11869a;, < x < 0.25 

^(a;,0) = 220, m = I 0.11869(0.5 - x), 0.25 < x < 0.75 ,v'^ = u^ = (22) 

[ 0.11869(x- 1), 0.75 < X < 1 

with the time varying distributed mechanical loading defined as 

F = 7000 sinMyj ^/ (mAm^) , G = 0. (23) 
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Under the given distributed mechanical loading, the SMA rod is expected to switch between 
different combinations of the martensite variants, and a hysteresis loop must be observed similar 
to those reported for ferroelastic materials at low temperature. In Fig. 2 we present simulation 
results for this case. The mechanical hysteresis is obtained by plotting displacement m as a 
function of F at a; = 3/8cm (the upper right plot). The time- varying mechanical loading for 
this case is plotted in the upper left plot. The simulated strain and the displacement distribution 
are also plotted as functions of time and space (lower plots). The combination of martensitic 
variants is changing with time-dependent mechanical loading and no stable austenite is observed 
at this low temperature. 

Our next goal is to analyze the behavior of the same SMA rod under a medium temperature 
where both martensite and austenite phases may co-exist. The following initial conditions will 
allow us to start from the austenitic phase: 

9{x, 0) = 250, uo = 0, v^^u^ ^ 0. (24) 

The boundary conditions as well as mechanical and thermal loadings in this case are kept 
identical the previous experiment. In this case, the free energy function has three minima 
that correspond to two martensites and one austenite. The numerical results for this case are 
presented in the left column of Fig. 3. It is observed that when the applied loading exceeds a 
certain value, the austenite is transformed to a combination of martensitic variants. The reverse 
transformation is taken place when the loading changes its sign. In contrast to the results 
presented in Fig. 2, we observe that the wide hysteresis loop, typical for the low temperature 
case, disappears. 

If we increase the initial temperature further to 6{x,0) = 300, the free energy function 
becomes convex and has only one minimum associated with the austenite phase. During the 
entire loading cycle, no martensite is expected under these thermal conditions. The dynamics 
of the SMA rod in this case exhibits nonlinear thermomechanical behaviour without phase 
transformations. This is confirmed by the numerical results presented in the right column of 
Fig.3. 



4.2 Thermally induced Phase Transformations and Hysteresis 

Thermally induced martensitic phase transformations and thermal hysteresis in SMA rods can 
be analyzed with the same model under time-dependent thermal loading conditions. Indeed, 
let us choose the initial conditions as follows: 

„, <,on / 0.11869a;, < x < 0.5 o i n 

e{x, 0) = 230, no = I Q^^ggg^^ _ < a; < 0.5 ' ^ = ^ = 

The boundary conditions remain the same as in the previous computational experiment, 
but the loadings conditions now become: 

G = 600 sin 9 1 {ms^cm) , F = 500 g/ {ms^cm^) . 

Numerical results for this case are presented in Fig.4. Analyzing strain and displacement 

distributions, we observe that the combination of martensitic variants is transformed into the 
austenite phase when the temperature exceeds a certain value. The reverse process is taken 
place when the temperature decreases, passing the critical threshold. Note that due to the 
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presence of thermal hysteresis, the critical temperature value for the martensite-to-austenite 
transformation is different from that of the austenite-to-martensite transformation. A schematic 
representation of the observed thermal hysteresis is given in the lower right part of Fig. 4 where 
we presented the temperature at a; = 3/8 as a function of strain at the same spatial point. 

5 Dynamics of SMA Patches 

The situation becomes more involved for 2D structures. Experimental, let alone numerical, 
results for this situation are scarce |2S]. In order to apply the FVM to the 2D model discussed 
in Section 2, we chose the same material as before, assuming that 02 = fci, as = ^2, 04 = /cs, ai = 
ki, as = 2ki and therefore effectively linking parameters in model and model (IE])- 

5.1 Nonlinear Thermomechanical Behavior 

The first numerical experiment on a SMA patch is aimed at the analysis of the dynamical 
thermo-mechanical response of the patch to a varying distributed mechanical loading, too small 
to induce any phase transformations. The initial temperature of the patch is set to 250° while 
all other variables are set initially to zero. Conditions at the boundaries are: 

0, Ml = 0, on left and right boundaries, 

(25) 

0, U2 = 0, on top and bottom boundaries. 

Similarly, the mechanical boundary conditions are enforced in terms of velocity components. 
The loading conditions in this experiment are: 

/i = 200 sin(7rt/6)g/(ms2cm2), /a = 200 sin(7rt/40)g/(ms2cm2). 

The time span for this simulation, [0,24], covers two periods of loading. The time stepsize is 
set to 1 X 10~^. We take 15 nodes used in each direction. The dimensions of the SMA patch 
are taken as 1 x 0.4 cm^. 

The variations in the displacements ui, U2 , deviatoric strain 62, and the temperature 6 
along the line y = 0.2cm (the central horizontal line) as functions of time are presented in 
Fig. 5. These simulations show that both thermal and mechanical fields are driven periodically 
by the distributed mechanical loading. Under such a small loading, the SMA patch behaves 
just like a conventional thermoelastic material. Observed oscillations are due to nonlinear 
thermomechanical coupling, but no phase transformations are observed in this case. 

5.2 Phase Transformations in SMA Patches 

Our aim in this section is to analyze spatio-temporal patterns of martensitic transformations 
in a 2D SMA patch. The SMA patch, used in this computational experiment, is made of the 
same material as before. The patch is assumed square in shape with dimensions 1x1 cm^. 
The initial temperature distribution is set to 6^ = 240°, and all other variables are initially set 
to zero. The boundary conditions are homogeneous and 

— = 0, U2 = Ui = 0, on all the four boundaries, (26) 



06 duo 

— = 0, — - 

dx ' dx 

89 ^ dui 

dy ' dy 
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where n is the unit normal vector. We apply the following loading to the sample, specified 
below for one period: 

sin(7rt/3), < t < 4, 

0, 10 < t < 12. 

The numerical results for this case are presented in Fig. 6 where the values of "y" are taken in 
the middle of the sample. The spatio-temporal plot of the order parameter 62 demonstrates a 
periodicity pattern in the observed phase transformations due to periodicity of the loading. It is 
observed also that the temperature oscillates synchronously with the mechanical field variables 
due to the thermo-mechanical coupling. 

As we mentioned earlier, there are two martensitic variants in the square-to-rectangular 
transformations. The following analysis proves to be useful in validating the results of compu- 
tational experiments. Assuming the temperature difference d6 = 6 — 60, one can easily calculate 
the deviatoric strain that corresponds to the austenite and martensite variants by minimizing 
the Landau free energy functional. In particular from the condition dFi/de2 = we get: 

r, 04 ± a/o-a — 4:a2d9aQ 

62 = 0; 62 = ^— . 

2afi 



The value 62 = corresponds to the austenitic phase. If we denote (04 + ^/a^'^^4a2d6a^) /2aQ 
by Cm, then 62+ = +^/e^ or 62- = —\fe^ are the strains that correspond to the two martensite 
variants. We call them martensite plus and martensite minus, respectively. If we take dQ = 42°, 
then for the material considered here we can estimate that 62+ = 0.12 and 62- = 0.12. This 
provides a fairly good estimate for the ID case. However, as was pointed out in [T2|, [15], for 
the 2D case such an estimate can be adequate only in homogenous cases. Although the quality 
of this estimate is dependent on the boundary conditions for a specific problem, this estimate 
proves to be a reasonable initial approximation to the deviatoric strain. 

In Fig. 7 we present two snapshots (at t = 2 and t = 8) of the spatial distributions of 62 
and 6. It is observed that when the mechanical loading achieves its (positive) maximum, the 
SMA patch is divided into two sub-domains determined by the deviatoric strain, as seen from 
the 62 plot at t = 2. In the upper- left triangular-shape area, the simulated deviatoric strain 
corresponds to the martensite plus, while on the opposite side, the deviatoric strain corresponds 
to the martensite minus. At t = 8, when the mechanical loading changes its sign to the opposite, 
the martensitic transformation is observed again, but now in the reverse direction. The second 
period of loading confirms these observations. 



6 Conclusion 

In this paper, we developed a finite volume methodology for the analysis of nonlinear cou- 
pled thermomechanical problems, focusing on the dynamics of SMA rods and patches. Both 
mechanically and thermally induced phase transformations, as well as hysteresis effects, in one- 
dimensional structures are successfully simulated. While these results can be obtained with 
the recently developed conservative difference schemes, their generalization to higher dimen- 
sional cases is not trivial. In this paper, we also highlighted the apphcation of the developed 
FVM to the 2D problems focusing on square-to-rectangular transformations in SMA materials 
demonstrating practical capabilities of the developed methodology. 
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Figure 1: Staggered grid for space discretization using finite volume methods. 




Figure 2: Mechanical induced phase transformation and mechanical hysteresis in SMA rod. 
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Figure 3: Dynamics and mechanical hysteresis of SMA at medium and high temperature 
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Figure 4: Thermally induced phase transformation and thermal hysteresis in SMA rod. 
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Figure 6: Martcnsitic transformation in a square SMA patch caused by mechanical loadings in 
the X and y direction. 
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